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I. INTRODUCTION 


Many experimental techniques have been employed to investigate 
ionospheric properties; however, little attention appears to have been 
given to the use of ion acoustic -wave phenomena. Although naturally 
occurring ion waves have been detected in the ionosphere [ l] , the author 
is unaware of any attempts to purposely create and exploit ion-wave phe- 
nomena. This dearth of interest in artificially generated ion waves is 
undoubtedly due to the difficulty encountered in laboratory attempts to 
generate and detect ion waves. However, sufficient success in the gener- 
ation and detection of ion waves in the laboratory has now been attained 
[2, 3] to warrant the consideration of satellite-borne ion-wave experi- 
ments. The purpose of this note is to discuss some of the problems and 
possibilities pertaining to the use of ion acoustic-wave phenomena for 
satellite experimentation in the ionosphere. 

II. POTENTIAL APPLICATIONS 

Ion-wave measurements have several potential applications. An 
accurate measurement of ion-wave velocities can yield direct information 
about electron temperatures. Furthermore, measurements of ion waves 
can provide information about the direction of magnetic field lines and 
information about the relative motion between a satellite and the ionospheric 
plasma. A knowledge of electron temperatures is important not only in 
itself but for the determination of electron number densities. 

A frequently used technique for measuring electron number densities 
is the electrostatic plasma probe. The current, I e , drawn by an electro- 
static probe is given by the relation: 


where j is the average electron current density (amp/cm 2 ) in the plasma, 
A s is the area of the plasma sheath, and F is the fraction of the current 
crossing the sheath which intercepts the probe. The area of the plasma 
sheath can be determined by self-consistent analyses of Poisson's 
equation, while F can be calculated by assuming that the plasma sheath 
edge and plasma probe form a diode [4]. The term j is given by the 
relationship: 


/ kT \ 2 

je = n el4^J *2.5X10"n.(T.) 


' 5 " e < T / 


( 2 ) 


where n e is the electron number density and T g is the electron temper- 
ature. If T can be determined independently and j has been measured 
with the probe (equation 1), then equation (2) can be employed to calcu- 
late the electron number density directly. 

The speed of propagation of an ion acoustic wave through a 
plasma may be specified as: 


(\^e\ 

v-hsrJ 


(3) 


where v^ is magnitude of the ion-wave velocity, Y e is the ratio of specific 
heats, k is Boltzmann's constant, and m^ is the ion mass. The atomic 
constituents of the ionosphere are well known and are displayed in Fig- 
ure 1. From Figure 1 it is seen that atomic oxygen is the predominant 
element above the altitude of 400 km. Consequently, the magnitude of 
the ion mass can be determined and, if the ion-wave velocity can be 
measured accurately, the electron temperature can be calculated. 

The application of ion- wave velocity measurements to the deter- 
mination of the relative motion between a satellite and the ionospheric 
plasma requires that the ion-wave velocity w be of the same order of 
magnitude as the satellite velocity v . Moreover, in order that the 
direction of the velocity vector may be determined, it is necessary to 
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make simultaneous measurements of ion acoustic velocity both in the 
direction of the satellite's orbital path and perpendicular to the satel- 
lite's orbital path. A quadrupole array of two perpendicular sets of 
transmitting -receiving plates which would permit orthogonal measure- 
ment of ion-wave velocities is shown in Figure 2. To measure the 
satellite velocity v g by the method illustrated, the array must be rotated 
until the planes of plates a and a' are parallel with the satellite's velocity 
vector. If v p is defined as the ion-wave velocity parallel to the orbital 
path, and v n a.s the ion-wave normal to the orbital path, then the satel- 
lite velocity is defined as; 



An experimental arrangement of equipment patterned after the 
quadrupole array shown in Figure 2 is discussed in greater detail later 
(page 9). However, it should be noted at this point that it will normally 
be necessary to rotate these electrodes in order to determine the direc- 
tion of the orbital velocity vector. 

The formula for the magnitude of the ion-wave velocity given 
in equation (3) is derived under the assumption that the electron tem- 
perature is much greater than the ion temperature. Because of the 
importance of this formula, a brief derivation is given in the section 
following. 


III. ION ACOUSTIC WAVE VELOCITY 

The correct expression for the ion acoustic -wave velocity has been 
the subject of some controversy, with most of the uncertainty centered 
upon the correct value of y . A simple derivation of an ion-wave equa- 
tion is given by Spitzer [5]; it begins with the equations of motion for 
ions and electrons: 

9 v i 9 v e 

Pi — = - VPi . and p e — = -vp e . (5) 

3 






If we assume an adiabatic pressure variation, the electron pressure may- 
be expressed as: 


p = Y k T n 

tr & 'o p « 


( 6 ) 


e 'e " e"e 

and there is a corresponding relation for the ion pressure. A correspon* 
dence between number density and velocity may be established through 
the equation of continuity: 


dn 

( 

Tt 


m* 


(pjv ) = - — (p Vv + v • Vp ) , 
'e e' e e e e 


where Z represents the degree of ionization of an atom. For small 
oscillations, the quantity v Q • Vp e may be neglected since it is the pro- 
duct of two first order terms. Hence, the equation of continuity for 
electrons reduces to: 

9^ 


a t 


m- 


pvt , 


( 7 ) 


and there is a corresponding relation for ions. If we differentiate equa- 
tion (6) and combine the result with equation (7), we obtain the general 
result for electrons or ions: 


at 


(1 + Z) Yk_T - v . „ 


m. 

x 


where it is assumed that the ion and electron kinetic temperatures are 
equal. Finally, if we combine this expression with equations (5) and 
derive a more accurate average over electrons and ions we obtain: 


92 v x = / Z ?e kT e + V T i \ 3v 
a t 2 l m, / a x 2 


x 


( 8 ) 


Equation (8) has the well known form of wave equation satisfied by a 
wave traveling in the x direction with a velocity obtained from the 
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equation: 


Z Y e k T e + Yj k Tj 

V C = * , 

m i 

It is shown below that unless T>> T, , ion waves are strongly damped; 
accordingly, the above expression for the ion wave velocity reduces to 
that given by equation (3). 

The magnitude of y has not been well defined since it is difficult 
to ascertain whether ion vibrations are isothermal or adiabatic. Essen- 
tially, y depends upon the number of degrees of freedom which the 
plasma particles may have, and is given by the relation: y = (2 + m)/m. 
Spitzer [5] argues that m = 1 in a fully ionized, compressionless plasma, 
and thus y should be 3. Stix [6] contends that ion oscillations are 
adiabatic (m = 3), and thus 7 should be 5/3. Experimental investiga- 
tions [3] appear to indicate that 7 should actually be closer to unity; 
hence, 7 will be assumed to have unity value for the calculations which 
follow. 

IV. LANDAU DAMPING OF ION ACOUSTIC WAVES, 

The probability of detecting ion waves depends to a large extent 
upon the amount of damping present. Unfortunately, theoretical con- 
sideration of collisionless damping is quite complicated since it involves 
the analysis of dispersion relations in the complex frequency plane. Be- 
cause of its importance, and for completeness, a brief discussion of the 
theory of collisionless damping of ion waves is given here. 

In 1945 , Landau [ 7 ], through an analysis of the linearized, 
collisionless Boltzmann equation, discovered that under certain condi- 
tions a ,! collisionless* ' damping of plasma oscillations could occur. Sub- 
sequent to Landaus work, a number of investigators studied the problem 
of Collisionless'' damping. One of the best treatments of ion acoustic 
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waves is contained in a paper by Fried and Gould [8], Fried and Gould 
solve the collisionless Boltzmann equation along with the Poisson equa- 
tion to obtain a dispersion relation of the form: 

x - 1 tz'd) + ez'dVe/?)] = o, (9) 

where 

£ ~ ~~ = ~~ [Re(to ) ± ilm(u) ) ] = x + iy; Q = T/T.; w is 

a a ^ X 

the complex frequency; k is the wave number; a is the mean thermal 

velocity; T e and T. are the electron and ion temperatures, respectively; 

■* • 

and Z ( § ) is the derivative of the "plasma dispersion relation. " A dis- 
cussion of the derivation of equation (9) is given in the appendix. Equa- 
tion (9) can be mimerically solved for £ by an iterative process in 
which values of x and y are chosen such that equation (9) is satisfied. 
Figure 3 exhibits solutions of equation (9) for 6 = 1,4, and 2,5. It can 
be seen that as 6 increases, the ratio Re(w)/Im(w) decreases, implying 
that for v lues of 0=3, 4, or greater, ion waves should be observable. 

In the laboratory experiments performed by Alexeff and Jones [3], 
electron temperatures as high as 10 °K were obtained, hence the require- 
ment for T e » was satisfied. Unfortunately, in the ionosphere T e is 
only about twice T ; hence, ion waves probably would be heavily damped. 
However, experimental results with cesium plasmas [z] in which T 0 and 
!• are approximately equal indicate that even for temall ratios of 0 , ion 
acoustic waves still can be detected. In this case, however, it is not 
certain how applicable equation (3) may be. 

Fried and Gould show that, if the ions and electrons are drifting 
relative to each other, ion acoustic oscillations may be enhanced. To an 
observer on a satellite, the ions have a undirectional drift velocity with 
respect to the satellite, while the electrons are essentially isotropic. 
Consequently, the motion of a satellite through the ionosphere may en- 
hance the detection of ion waves. 
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V. SATELLITE ION WAVE EXPERIMENTS 


Satellite -borne ion-wave experiments will have to be planned 
carefully since there may be complicating factors. The ionospheric 
plasma differs considerably from a quiescent laboratory plasma in that 
the Earth's magnetic field and perturbations caused by the passage of 
the satellite will affect the behavior of the ion waves. Furthermore, when 
ion waves are generated, perhaps by a voltage pulse, it is possible to 
excite electromagnetic waves which may tend to confuse experimental 
data. 


A. Theoretical Considerations 

The probable ion-wave velocity spectrum of the ionosphere can 
be established after the electron temperature range has been determined. 
Figures 4 and 5 show the approximate electron temperatures and electron 
concentrations in the ionosphere for altitudes 100 km to above 700 km [9] . 
The acoustic velocities for waves in the ionosphere may be determined 
by substituting the appropriate values obtained from Figure 4 into equation 
(3), For proton ion waves (m^ = 1 amu), the velocity is expressed as: 


v s 9 IVt" meter /sec, 


( 10 ) 


while for any other atom the velocity is expressed as: v = v^VA , where 
A is the atomic number. Shown in Figure 6 is a range of acoustic veloci- 
ties for protons, oxygen, and nitrogen. From Figures 4 and 6 we can 
deduce that at an altitude of approximately 500 km, where T e ^2 X 10 3 °K, 
an oxygen acoustic wave would have a velocity of 1. 0 km./ sec. Typical 
satellite velocities are in the range of 7 to 8 km/ sec; consequently, the 
ion wave velocities have a velocity sufficiently comparable to the satellite 
velocity to be detectable. 

For situations in which dispersion effects can be ignored, the 
angular frequency of an ion wave is given by the equation: 


2 Zm p 2 

03 ? ~ ® 03 

1 m; P ' 


( 11 ) 
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in which 


_ 


4irne z 

c 

, and co is the electron plasma frequency, 

m e P 


In general, the ions encountered in the ionosphere are singly 
ionized; that is, Z = 1. The numerical value of f p = W / 2 7T is 
8, 97 X 10 3 (nj 2 ) [5] ; therefore: 


f. = 209(n G /A) 2 hertz (Hz). 

X ^ 


(12) 


Thus, the ion wave frequency for protons is: f = 209 ( n j) Hz, for 
nitrogen it is: f = 56 (n^) Hz, and for oxygen it is: f = 52. 2 (n^J Hz. 
Ion wave frequencies corresponding to the electron densities depicted 
in Figure 5 are shown in Figure 7, 

The Earth’s magnetic field exerts an influence upon the ion waves 
and, therefore, may affect their detectability. The parameters which 
determine the magnitude of the influence are the ion cyclotron frequency: 


co B Z£ B 
C m . 


(13) 


and the ion gyro radius: 




(14) 


in which 33 is the magnetic field intensity in gauss and is the component 

of perpendicular to B_. The cyclotron frequency f = <o /2ir has the 

o c 

numerical form f = (1. 54X 10 3 ZB)/A, The Earth's magnetic field at 
satellite altitudes is approximately 0.5 gauss [ll] ; thus, for singly ion- 
ized atoms the cyclotron frequency is; 


f ~ 0. 77A” 2 kilohertz (kHz), 
c 


(15) 


Hence, i Q ~ 0. 77 kHz for protons* 0. 206 kHz for nitrogen ions, and 
0. 193 kHz for oxygen ions. As indicated by Figure 7, the acoustic fre- 
quency will be much higher than the ion cyclotron frequency; therefore, 
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there should be little interference between the cyclotron motions and 
the acoustic waves from a frequency standpoint. 

The cyclotron radius has the numerical form: 

Ao> x 

r = 1. 04 X 10~ 2 — - — centimeters. (16) 

Jd 

In a 0. 5 gauss magnetic field, r = 16 cm for protons, 60 cm for nitrogen, 
and 64 cm for oxygen. In order to assess the influence of the ion cyclotron 
radius, we must ascertain the wavelength of the ion acoustic waves, 

X = v/f. The wavelength may be obtained from equations (10) and (11) in 
the form: 

JL 

X = 43. 5(T e /n e ) 2 centimeters (17) 

for which X is independent of A. A graph of wavelength versus altitude, 
in which data from Figures 4 and 5 are used, is given in Figure 8. From 
Figure 8 we see that at 500 km the ion acoustic wavelength is 4. 35 cm. 
From this discussion of wavelengths we conclude two things. 

First, the acoustic wavelength is of a magnitude commensurate with 
reasonable spacecraft dimensions and, second, the acoustic wavelength 
is smaller than the ion cyclotron radius. Consequently, the ion waves 
will appear as a perturbation on the ion cyclotron motion and the, cyclotron 
motion should not limit the detectability of the ion waves. 

B. Experimental Apparatus 

The experimental apparatus for creating and detecting ion waves 
can be patterned after the equipment employed by Alexeff and Jones [3] . 
Shown in Figure 9 is a possible experimental arrangement of the quadru- 
pole transmitter-receiver mentioned on page 3 to create and detect waves 
perpendicular to each other. The apparatus should be constructed to per- 
mit variable plate separation, and the entire receiver should be able to 
rotate with respect to the satellite. The latter capability is important 
when measurements are being made of the satellite velocity vector relative 
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to the plasma. Variable plate separation is important because the damping 
factor is uncertain, Wong et al [2] show that when and are approxi- 
mately equal, an acoustic wave is effectively damped in approximately 1/4 
wavelength ( 1 / e) ; however, for best wave resolution, plate separation 
should be at least one wavelength. 

An ion wave having a velocity of 1 km/ sec traverses 0, 1 cm per 
fisec. Thus for a plate separation of one wavelength (^5 cm), it takes 
an ion wave approximately 50 fJ, sec to pass between the plates. A concep- 
tion of a time - of-flight circuit capable of measuring ion waves with this 
magnitude of transit time is illustrated in Figure 9. Basic timing is pro- 
vided by a crystal-controlled oscillator. The output of the oscillator is 
differentiated to provide excitation pulses at the transmitting plates and 
for the "one- shot" multivibrators OS-1 and OS-2. Both "one- shots" have 
varible pulse widths, which permit the direct measurement of the transit 
time of a wave between plates. The pulse from OS- 2 acts as a timing gate. 
If the received pulse and the OS-2 pulse are coincident, .the "and" circuit 
will respond and the received pulse will be detected. Varying the pulse 
width of OS-1 will give the OS-2 pulse the proper position, in time, to be 
coincident with the received pulse. Initially, the OS-2 pulse should be 
fairly wide; then, when the received pulse is centered within the OS- 2 
pulse, the OS-2 pulse can be narrowed to refine the measurement of the 
acoustic-wave transit time. Resolution on the order of 1 /Xsec should be 
attainable over a range of approximately 1 msec. Thus, a velocity incre- 
ment resolution of about 1 cm/sec should be possible. 
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FIGURE 2. RELATIONSHIP BETWEEN SATELLITE VELOCITY 
AND ION ACOUSTIC-WAVE VELOCITY 
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FIGURE 3. SOLUTIONS OF EQUATION 9 ILLUSTRATING THE 
EFFECT OF 6 = T^/T. UPON ION ACOUSTIC- WAVE DAMPING 
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FIGURE 5. DAYTIME IONOSPHERIC ELECTRON CONCENTRATIONS 
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FIGURE 6. ION ACOUSTIC- WAVE SPEEDS FOR TYPICAL 
IONOSPHERIC ELECTRON TEMPERATURES 
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FIGURE 7. ION ACOUSTIC- WAVE FREQUENCIES FOR 
TYPICAL IONOSPHERIC ELECTRON CONCENTRATIONS 
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APPENDIX 


DERIVATION OF LONGITUDINAL OSCILLATION 
DAMPING FACTOR 

One of the more rigorous approaches to the analysis of plasma 
oscillation phenomena involves the use of the plasma dispersion rela- 
tion applicable to the plasma situation under investigation. Plasma 
dispersion relations have been investigated by many authors; the following 
analysis is a recapitulation of the work of Fried and Gould [8] and 
Jackson [10 ] , 

We begin with the collisionless Boltzmann equation for electrons: 


— + v . V^f - — (E + vXB)vJ = 0 (la) ' 

St r m v 

v e 

where f(:r, v, t) is the electron distribution function, r’ and v are the 
displacement vectors in configuration and velocity space, e/m is the 

V/ 

ratio of electron charge to mass, V is the gradient operator, and E and 
— * 

B are the electric field and magnetic field intensities. Since we are 
interested in. longitudinal oscillations, the cross-product of v and B is 
assumed equal to zero. We shall further assume that the oscillations 
are one -dimensional and that the velocity space distribution function is 
spherically symmetric; cherefore, equation (la ) has the form.: 


d± 0f 

at v a t 


m. 


E |i = 0 . 

3 V 


Next let us assume that f and E have the form: 


f(x, v, t) - f 0 (v) + fj (x, v, t) 


(2a) 


E(x, v, t) = E 0 (v) + Ej(x, v, t) 

with f 1 «f 0 and E 1 «E Q , that is, fj and E i are small disturbances to the 
■ equilibrium quantities f Q and E Q . Substituting into equation (2a) we 
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obtain: 


di. e / 01 at \ 

- + + 37) =° 

since by definition df^/dt - 3f 0 /3x = 0. The electric field E Q can be 
omitted with no loss of generality since E Q represents an unchanging 
field and we are interested only in the behavior of the perturbed field 
E, In addition, we shall neglect second order (nonlinear) effects, that 
is, E t 3^ 3v. Thus we arrive at the linearized Boltzmann equation for 


electrons: 


0 f l . 0 f ! 

sT +v “37 


= 0 . 


(3a) 


A similiar equation can be obtained for ions, with F replacing f and 
+ eZ/m. replacing -e/m e< We shall assume that Z = 1,0 in the following 
analysis. 

The coupling between the electric field and the distribution 
function arises through the Poisson equation, which for one dimension 


has the form: 


J OO 

- = 4 7Tp = -e47T (n* - n ) = 4tt e f dv(f - F ), 

1 c -00 1 1 


(4a) 


where p is the electric charge density. 

Equations (3a) and (4a) are conveniently simplified by employing 
the Fourier and Laplace transforms defined by: 


and 


F [f(t)] = J exp(ico t)f(t)dt = h( w ) 


L [f(t)] = J exp - (pt)f(t)dt = g(p) 


(5a) 
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with the inversion relations 


a+ioo 


“l 1 c 

h [ g<p) 3 = 2^^ J . ex P(pt)g(p)dp 


a-ioo 


(6a) 


F [h(co)]= J exp - (ico t)h(co )dco . 


Applying the Fourier transform to x (integrating by parts and observing 
that ^(±00= 0)), we obtain: 

r°° 9f i 

J v exp(ikx)dx 


00 w 

v l f i exp(ikx) - ivk J £^ exp(ikx)dx 


= -ivkf x (k, v, t) . 

Applying the Laplace transform to t (setting p = -fico), we obtain: 

r° 3f i 

J -0 exp -(io)t)dt = icof 1 (k, v, o ) 
equation ( 3 a) then becomes: 


io,>f,(k, v, o>) - ivk • f x (k, v, 00) « 


J0v 


from which we obtain 


f, = — 

1 m. 


e i“ 57 (i " ‘ ivk) " • 


(7a) 


In an analogous manner' we obtain Fj : 


p = E (ico - ivk) 

1 m • lev 


Finally, we Fourier transform Poisson's equation to obtain: 


-ikE = 4 ?re f dv(F -f ) . 

1 V 11 


(8a) 


Substituting f and F into equation (9a), we obtain: 


kEj = ^ i£) (v » r i , 

1 km e 4o 1 \ 0v m. d v / ' k 7 


m. 

1 


which is rearranged to: 


E 


.[ 


k 2 _ 4*e_ na. 


m. 


00 / , 

/ T 

-00 \ 0 


1 3£( 


m, 


+ 


m i n o 


£5 

a v 




where n is the particle number density. Some simplification is possible 
by noting that the "plasma" frequency relation is: 

47rne 2 /m = 

C JJ 

and by introducing: 

6 = m e /m i , 

whereupon we have the dispersion relation: 

r °° / 91 3F \ 

-/ k + 6 tt) < v 


O'/I to \ — ^ 

H (k, 'k> = -C2 


-00 


U) 


-1 


) = 0 


(10a) 


In principle, if f 0 and are known, equation (10a) can be 
solved for a) , and it is seen that _co_ is a function of velocity, that is, 
the frequencies are dispersed in velocity space. In general, co is a 
complex number and it is well known that if Im(to)< 0, a wave of the 
form exp - (icot) will be damped, while if Im(oo)> 0, a wave will grow and 
unstable oscillations can occur. The question of collisionless damping, 
usually known as Landau damping, has been debated considerably. A 
principal difficulty arises because, as Jackson demonstrates [10], the 
dispersion relation of equation (10a) is valid only for Im(o))> 0, and to 
obtain a relation for Im( co ) < 0, equation (10a) must be analytically con- 
tinued across the real axis. 
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In the case of Maxwellian distribution?, analytic continuation 
is possible; therefore, we shall consider these types of distributions. 
Accordingly, we assume f Q and have the form: 


f (v) = eX P- (V i 2/a2) 


and 


a7r‘ 


F (v) = 6X P- (v2 / A2) 
A 7 r 2 


(11a) 


where the relation a 2 = 2kT e /m is the mean thermal speed squared. 
Substituting equation (11a) into equation (10a), we have 


0 = 


“p 


00 r 

+ / dv [ 

- oo 


2v 

a 


r* exp 


2v 

+ 6 T 




a tt ‘ 


exp 


- 6 ) 


•(5)i (’--a 


-l 


This expression may be reduced if we employ: 


03 


m 2 
e _a 

m- a 2 


T i 


, and x = v/a . 


Therefore: 


0 = / dx[xexp-(x 2 ) 

P a7T2 -00 


+ 0 X exp - (x 2 0 /5 ) ] |xa - 


- 1 


The term w is a well defined analytic function and may be ex 
tended into the complex plane when the following definition is used: 


£ = x h iy = u>/ka . 
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Thus the dispersion relation is reduced to: 


OO 


0 = 


to 2 


a Z TT Z -oo 


/ dx [x exp - (x 2 ) 


"* 1 

+ 0 x y/6f 6 exp - (x 2 0 / 6) ] [x - £ ] . (12a) 


The Debye wave number is given as follows: 
k D = 2 “p /a2 


and if we define 


we obtain 


X = k 2 /k^ 


oo 

0 - X + 7r 2 J" dxx [ exp - (x 2 ) 


— oO 


+ 0V0/6 exp - (x 2 9/ 6)] [x - £] 


- 1 


(13a) 


It is now convenient to consider the plasma dispersion function 
-1 00 


Z(z) = 7r 2 f du exp - (u 2 ) (u - z) 1 , 


(14a) 


-00 


Equation (13a) is in the form of a Cauchy integral J"dzf(z)/(z - z Q ), and 
from the theory of complex variables the derivative of this integral 
can be written as; 


d r f(z)dz _ r f'(z) 
dz ^ z - z n z - z„ 


(15a) 


-0 ~ “0 

(prime denotes differentiation). Thus the relation dZ(z)/dz = Z (z) 
becomes: 

, - k r° 

Z (z) = -2 7 r 2 J du u exp - (u 2 ) (u - z) 

-00 

When equation (15a) is compared with equation (13a), it is seen that 
equation (13a) reduces to: 

x- §[z'< i) + e z'dVeTT)] = o. < 16a > 
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Although equation (16a) may be solved for £ , the process is 
laborious at best. Fortunately the function Z(£) closely resembles the 
error function for complex argument defined by: 

W(z) s — r f du exp - (u 2 ) (u - z) 1 

7TI J 
-oo 

From a comparison of the definition of W(z) with equation (14a), we 
obtain: 

z(£) =v^ w(g)/i . 

The derivative of W(z) is given by the relation: 

W'(z) = -2zW(z) + 2i7r” 2 . 

Since the dispersion relation is symmetric in x, we can obtain values 
for £ = x-iy by employing the relation: 

W(-z) = 2 exp - (z 2 ) - W(z) . 

Employing these relationships, we find: 

z'(£) = [- £WU)+ — U- j 

it vf J 

= 2 [1 + i Vir £W(g)] 

and 

z'(-6) = 2 [1 + iv9 e (2exp-(e 2 ) - W(g)JJ . ! (17a) 

If equation (17a) is substituted in equation (16a) and expanded, 
numerical results can be obtained for g = ook/a. The procedure, while 
tedious, is straightforward. One first determines those values of x 
and y which reduce the imaginary part of equation (16a) to zero. Then 
the real part of equation (16a) determines X . This procedure has 
been carried out by Fried and Gould, and the resultant curves in the £ 
plane are exhibited in Figure 3. It should be noted that as 6 increases, 
Im(oo) becomes quite small, indicating that 0 2:2 or 3, there is a reason- 
able probability of observing longitudinal plasma waves. 


